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We revisit the quantum / classical integrable model correspondence in the context of 
inhomogeneous finite length XXZ spin- = chains with periodic boundary conditions 
and show that the Bethe scalar product of an arbitrary state and a Bethe eigenstate 
is a discrete KP r-function. The continuous Miwa variables of discrete KP are the 
rapidities of the arbitrary state. 



1. Introduction. 

Quantum models of the statistical mechanical type (the only quantum 
models discussed in this note) such as integrable 1-dimensional quantum 
spin chains, and classical models such as integrable nonlinear partial differ- 
ential equations, are related in the sense that the methods used to compute 
in the former, particularly the quantum inverse scattering transform, also 
known as the algebraic Bethe Ansatz, are quantum versions of those used 
to compute in the latter, namely the classical inverse scattering transform. 
It is therefore natural to expect that the quantum integrable models have 
classical limits in which they reduce to classical counterparts. 

What is less than natural to expect, at least to our minds, is that basic 
objects in quantum integrable models, such as the correlation functions, 
turn out to have direct interpretations in terms of objects in classical inte- 
grable models, such as solutions of integrable nonlinear partial differential 
and difference equations, without taking a classical limit. But this turns 
out to be the case, and it points to a direct connection between quantum 
and classical integrable models that is distinct from, and to our minds at 
least as fundamental as that obtained by taking a classical limit. 
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Notes on the literature. The following is far from a comprehensive sur- 
vey of the relevant literature. To the best of our knowledge, a direct con- 
nection between quantum (statistical mechanical) and classical models of 
the type that we are interested in first appeared in [1], where Ising spin- 
spin correlation functions in the scaling limit were shown to satisfy Painleve 
equation of the third kind, and subsequently in [2,3], where critical Ising 
correlation functions on the lattice were shown to satisfy the Toda lattice 
equation in Hirota's bilinear form. Further results, along the same lines 
as in [1], for the XXZ spin chain at the free fermion point, were obtained 
in [4], as reviewed in [5]. 

The fact that r-functions (solutions of Hirota's bilinear equations) ap- 
pear in the Ising model as well as in KP theory was discussed in works by 
the Kyoto group and reviewed in [6] where it was argued that the mathe- 
matical reason underlying this coincidence is the fact that both quantum 
and classical models are based on infinite dimensional Lie algebras that are 
realized in terms of free fermions. 

Closest to the spirit of this note is the work of Krichevcr et al. [7] , re- 
viewed in [8]. The starting point of [7] is the observation that the Bethe 
eigenvalues satisfy Hirota's difference equation, various limits of which lead 
to a large number of integrable differential and difference equations [9]. 
We will comment on the results of [7] and how they differ from the re- 
sult in this note in section 6. More recently, studies of the ultra-discrete 
limit of quantum integrable spin chains revealed many classical integrable 
structures [10]. 

Finally, while we are only interested in integrable quantum models in 
statistical mechanics in this note, it is important to mention bosonisation 
(the operator formulation of Sato's theory) as a deep and established cor- 
respondence between the quantum field theories of free fermions, which are 
integrable quantum models, and classical integrable hierarchies, as reviewed 
in [11]. In this correspondence, expectation values of fermion operators have 
direct interpretations in terms of solutions of integrable nonlinear partial 
differential equations. Bosonization was further extended to connect KP 
theory and conformal field theories on Riemann surfaces (which are inte- 
grable quantum models) in [12]. 

The long term aim of our work is to develop a correspondence between 
integrable statistical mechanical models and classical integrable hierarchies 
that is as direct and detailed as that obtained by bosonisation between free 
fermions and classical integrable hierarchies. 
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Bethe scalar products and continuous KP r-functions. Consider 
the inhomogeneous length- L XXZ spin-^ chain with periodic boundary con- 
ditions. Following [13], the Bethe scalar product (Ai, • • • , \n\hi, ■ • • , Atjv)/3 
of an arbitrary state (Ai,--- , Ajv| where the auxiliary space rapidities 
{Ai, • • • , Aat} are free, and a Bethe eigenstate - - - , /an) f) where the aux- 
iliary space rapidities {/xi,-- - , /ijv} obey the Bethe equations, is a poly- 
nomial T-function of the continuous (differential) KP hierarchy. In this 
identification, the KP time variables {ti, • ■ ■ } are power sums of the free 
rapidities {Ai, • • • , Aat}. However, these polynomial KP T-functions involve 
by construction more time variables than free rapidities. The reason is as 
follows. 

Expanding the scalar product in terms of Schur polynomials s\, as- 
sociated to Young diagrams {A}, that are functions of the rapidities 
{Ai,-- - ,Aat}, the maximal number of rows in any Young diagram A is 
N . Switching to KP time variables {ti, t 2 , • • ■ } that are powers sums in the 
rapidities, we obtain character polynomials x\ that depend on effectively as 
many time variables as the number of cells in (that is, the size of) A which 
is larger than N. Consequently, the KP time variables {t\, t 2l - • ■} were for- 
mally considered in [13] to be independent, and the Bethe scalar product 
was defined as a restricted KP r-function obtained by setting {t\ , ti , ■ ■ ■ } to 
be power sums of a smaller number of independent variables {Ai, • • • ,Aat}. 

Bethe scalar products and discrete KP r-functions. In this note, we 
simplify the correspondence of [13] by working solely in terms of the free 
rapidities {Ai, • • • ,\n} which are now continuous Miwa variables and the 
T-functions that we obtain are those of the discrete KP hierarchy [14, 15]. 

Outline of contents. In section 2, we recall basic facts related to sym- 
metric functions, Casoratian matrices and Casoratian determinants. In 3, 
we recall basic facts related to the XXZ spin-i chain, the algebraic Bethe 
Ansatz, the Bethe scalar product, recall Slavnov's determinant expression 
of the Bethe scalar product and show that it is a Casoratian determinant. 
In 4, we recall basic facts related to the continuous and discrete KP hierar- 
chies and define the Miwa variables that relate the two. In 5, we show that 
Bethe scalar products in the XXZ spin-i chain with periodic boundary con- 
ditions are discrete KP r-functions. In 6, we collect a number of remarks. 
Space limitations allow us to give no more than the minimal definitions 
necessary to fix the notation and terminology supplemented by references 
to relevant sources. 
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2. Symmetric functions and Casoratians. 

The canonical reference to symmetric functions is [16]. Casoratian matrices 
and determinants are carefully discussed in [15]. The definitions in [15] are 
more general than those used in this note. 

Frequently used notation. We use {x} for the set of finitely many vari- 
ables {xi, x 2l ■ ■ ■ , %n}, or infinitely many variables {x\, X2, ■ ■ ■ }■ The car- 
dinality of the set should be clear from the context. We use {x m } for 
{x} but with the element x m missing. In the case of sets with a re- 
peated variable x i7 we use the superscript (m,) to indicate the multiplic- 
ity of Xi, as in x^ n '\ For example, {xf\ x 2 , x^\ x±, ■ ■ ■ } is the same as 
{xi,x\,Xi,X2,X3,X3,Xi, ■ ■ ■} and /{• • • , x[ mi \ • • • } is equivalent to saying 
that / depends on distinct variables all of which have the same value 
Xi. For simplicity, we use Xi to indicate x^p . In calculations, it is safer 
to think of any Xi with multiplicity rm > 1 initially as distinct, that is 
{xi,i, %i,2, ■ ■ ■ ,Xi. mi }, then set these variables equal to the same value 
Xi at the end. 

We use the bracket notation [x] — e x — e~ x , and 

A{x}= H (xi-xj), A t „ 9 {A}= J] (!) 

l<i<j<N l<i<j<N 

for the Vandermonde determinant and its trigonometric analogue. 

The elementary symmetric function ei{x} in TV variables {x} is the 
coefficient of k % in the expansion 

N oo 

Y[{i+xik) = Y, e ^ kl ( 2 ) 

i=l i=0 

For example, e {x} — 1, e 1 (xi,x 2 ,x 3 ) = X\ + x 2 + x 3 , e 2 (x 1 ,x 2 ) = X\X 2 - 
ei{x} = 0, for i < and for i > N. 

The complete symmetric function hi{x} in N variables {x} is the 
coefficient of k % in the expansion 

N oo 

I1t^ = £mzK (3) 

For example, h {x} = 1, hi(xi,x 2 ,x 3 ) = Xi + x 2 + x 3 , h 2 (x u x 2 ) = x\ + 
x\x 2 + x 2 , and hi{x} = for i < 0. 



May 24, 2010 1:33 Proceedings Trim Size: 9in x 6in 



Miwa-3 



Useful identities for hi{x}. From Equation (3), it is straightforward to 
show that 



hi{x} = hi{x m } + x m hi-i{x} (4) 
and from that, one obtains 

h i {x 1 ,x 2 , ■ ■ ■ ,x N } = 

hi{xf\x 2 , ■ ■ ■ ,x N } - x 1 h i - 1 {xf\x 2 , ■ ■ ■ ,x N } (5) 

(x 2 - Xi)hi{xf\x^\x3, ■ ■ ■ ,x N } = 

x 2 hi{xi,x^\ ■ ■ ■ ,x N } - xih i {x ( j 2 \x 2 , ■ ■ ■ ,x N } (6) 

The discrete derivative A. m hi{x} of hi{x] with respect to any one 
variable x m £ {x} is defined using Equation (4) as 

Amhi{x} = feiM-M^m} = hi i{x} (7) 

Xm 

Note that the effect of applying A m to hi{x} is a complete symmetric 
function hi-i{x} of degree i — 1 in the same set of variables {x}. The 
difference operator in Equation (7) is not the most general definition of a 
discrete derivative, but it is sufficient for the purposes of this note. For a 
more general definition, see [15]. 

The Schur polynomial sa{^} indexed by a Young diagram A = 
[Ai, . . . , A r ] with Xi ^ 0, for 1 < i < r, and Ai = 0, for r + 1 < i < N, is 

det (x X ^ 3+N ) 

s,{x } = 1 A{x > ^ = det (JW,)) ^ (8) 

For example, s<p{x} = 1, S[i](xi, x 2 , x 3 ) = xi+x 2 +x 3 , S[i,i](a;i, x 2 ) = X\x 2 . 
The first equality in Equation (8) is the definition of s\{x}. The second 
is the Jacobi-Trudi identity for s\{x}. s\{x} is symmetric in the elements 
of {x} and requires no more than r (the number of non-zero rows in A) 
variables to be non- vanishing. 
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The one-row character polynomial indexed by a one-row Young 

diagram of length i, is the i-th coefficient in the generating series 



5^Xi{t}** = exp 



(9) 



For example, xo{t} = 1, XiW = h, X2{t} = y +*2, = ^+*i*2+*3, 

and Xi{t} = for « < 0. Since U has degree i, \i is n °t symmetric in {t} 
and generally depends on as many t- variables as the row- length i. 



The character polynomial Xa{*} indexed by a Young diagram A = 
[Ai, A 2 , . . . , A r ] with r non-zero- length rows, r < N, is 



XA{*} = det (xxi-i+jlt}) (10) 

t 2 t 3 t 4 

For example X[i,i]{*} = 2"-*2, X[2,i]{*} = y-*3, X[2,2]{*} = t^-*i*3+*2- 
Notice that Xa{*} can depend on all t,, for i < |A|, where |A| is the sum of 
the lengths of all rows in (or area of) A. 

From character polynomials to Schur polynomials. Assuming that 
the ^-variables are independent and that we have sufficiently many x- 
variables, then setting t m — > ^ J2i=i X T sends Xi{t} ~^ hi{x}. In this 
note, as in [13], we study Bethe scalar products that are polynomials in TV 
variables {x\ 1 x%, ■ ■ ■ ,xn}. We can expand these scalar products in terms of 
Schur polynomials s\{x} where {A} has at most N rows, or in terms of the 
corresponding character polynomials X\{t} that require more t- variables 
(which are power sums in the x- variables) than TV and therefore cannot be 
all independent. We choose to work in terms of the x- variables and s\{x}. 

Casoratian matrices and determinants. A Casoratian matrix M of 
the type that appears in this note is such that the elements My satisfy 
either 



M iJ+1 {x} = A m My{x}, or M i+hj {x} = A m M i:j {x} (11) 

where the discrete derivative A m is taken with respect to any one variable 
x m G {x}. If M is a Casoratian matrix, then detM is a Casoratian de- 
terminant. Casoratian determinants are discrete analogues of Wronskian 
determinants. 
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3. The XXZ spin-i chain and the Algebraic Bethe Ansatz. 

The XXZ spin-i chain is discussed in detail in [17,18]. A standard reference 
to the algebraic Bethe Ansatz, including the Bethe scalar product and 
Slavnov's determinant expression, is [5] . We leave the definition of auxiliary 
and quantum spaces, auxiliary and quantum rapidities, and the precise 
action of the various operators to [5] . 

Frequently used variables. In the following, L is the number of sites in a 
periodic XXZ spin-i chain, and N is the number of Bethe operators B(/ii) 
that act on the reference state |0) to create an XXZ state , • • • A 
is also the rank of the matrix whose determinant is Slavnov's expression for 
the Bethe scalar product. We use the set {A} for the free auxiliary space 
rapidities, {/i} or more explicitly {/ip} for the auxiliary space rapidities 
that satisfy the Bethe equations, and {v} for the quantum space rapidities 
(the inhomogeneities) . A Bethe eigenstate state whose rapidities satisfy 
the Bethe equations is also denoted by a subscript (3, such as \\}p. 7 is 
the crossing parameter. We use the exponentiated variables {xi, iji, Zi, q} 
— {e Ai , e Mi , e" 4 , e 7 }, but still refer to the exponentiated variables {x,y, z} 
as rapidities rather than exponentiated rapidities for simplicity. 

The i-operator of the XXZ spin-^ chain is 



L a i{\,v) = 



' [A - v + 7] 
\\-v\ [7] 

[7] [A - v\ 

[A - v + 7] , 



(12) 



where a is an auxiliary space index and i is a quantum space index. 



The monodromy matrix of the inhomogeneous length-!/ XXZ spin-i 
chain is 



r -<M3w^).-S i -<*-' <13) 

where it is conventional to suppress the dependence on the inhomogeneous 
quantum space rapidities Vi in T a and its elements, and each of the operators 
A, B, C, and D acts in the tensor product V\ ® • • • ® Vl where Vi is a vector 
space isomorphic to C 2 . 
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The transfer matrix is the trace of the monodromy matrix over the 
auxiliary space, 

Tr a T a (A) = A(X)+D(X) (14) 

An arbitrary state is generated by the action of the B((i) operators 
on the reference state |0) = <g) L ( o ) ' 

\n)=B(n 1 )...B(n N )\0) (15) 

An arbitrary dual state (A| is generated by the action of the C(A) 
operators on the dual reference state (0| = ® L ^1 oj , 

(\\ = (0\C(X 1 )...C(X N ) (16) 
The scalar product of a state and a dual state is 

<A|a»> - <0| C(Ai) . . . C(X N ) Bim) . . . B(fi N )\0) (17) 
A Bethe eigenstate \/J,)p is an eigenstate of the transfer matrix, 

[a{\)+D(X)}\^ p =E(XMp (18) 

where E(X) is the corresponding Bethe eigenvalue. For a state \fi) to be a 
Bethe eigenstate, its auxiliary space rapidities must satisfy a set of Bethe 
equations. 

The Bethe equations that must be satisfied by the N auxiliary space 
rapidities of a state = B(fii) ■ ■ ■ B([i N )\0) in order to be a Bethe eigen- 
state, in the specific case of the inhomogeneous lcngth-L spin-i chain, are 

nt>-^+7] A [^-M J -7] = 1 (19) 

where {vi, • • • , vl}, are the quantum space rapidities, which are taken to be 
part of the parameters that specify the spin chain, rather than the definition 
of the Bethe state. 
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A Bethe scalar product is a scalar product of an arbitrary state (A| and 
a Bethe eigenstate \n)p, 



(X\fi)p = (0| C(Ai) . . . C(X N ) B(ji 1<p ) . . . B( W )|0) (20) 

Bethe scalar products as in Equation (20) play a central role in computing 
XXZ correlation functions [19], hence their importance. 

Slavnov's determinant expression. In [20], Slavnov obtained an ele- 
gant determinant expression for the Bethe scalar product, 



TT^ [A - u + t1 n l 

{Mfi)fj = b]N a{a}am n n ^ ^ ^ ^ detn ( 2i ) 



where the components of the NxN matrix ft are 



^ [K - m +7] 

[Mj - A,] \m - Xi + 7] y [Ai - vh] f = \ [Xi -fn+j] 

Slavnov's scalar product is the main object of interest in this note. We wish 
to show that it is a Casoratian determinant and that the latter satisfy the 
bilinear identities of a discrete KP hierarchy [15]. 

Re- writing Slavnov's determinant expression. In [13], it was found 
useful to rewrite Slavnov's determinant expression for the Bethe scalar 
product as follows. First, we change variables and work in terms of ex- 
ponentials of the original variables as follows 



{e 2X %e 2 ^,e 2 ^,e^}^{x u y t ,z u q} (23) 

but continue to call the exponentials {x, y, z} rapidities as that is simpler 
and should cause no confusion. Ignoring prcfactors that do not depend on 
{x}, it was shown in [13] that the relevant part of Slavnov's determinant 
expression can be re-written as 
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detQ 



and 



, det ft 



N+L-l 



-, where fly = J2 x * * l = _ E fj * ^ 

^ W fc=i z=i 

(24) 



Pij 



k m— 1 
L 



JV 



v m— 1 



AT 



II ~ ^"3 2 ) • e (£+Jv-;){-&9 2 }{-2<? 2 } 

(25) 



In Equation (25), ek{yj}{z} is the fc-th elementary symmetric polynomial 
in the set of variables {y} U {z} with the omission of yj. 

A Bethe scalar product is a Casoratian determinant. We wish to 
show that Slavnov's determinant expression is Casoratian in the free ra- 
pidities {x} of the general state. Expanding detfi, using the Cauchy-Binet 
identity, we obtain 



det O = det 



( N+L-l 

x i 1 

*;=i 

E 

l<fei<---<fciv<W+L-l 



det ^x^ 3 1 j det ^K/^jj 
]T det (,-,^ ') det 

0<A N <---<Ai<L-l 

E dct (^ +JV+1 " J_1 ) det (/t Ai+A r+i_ i;j ) 

0<A N <---<Ai<L-l 



(26) 



From the definition of Schur polynomials that uses the Jacobi-Trudi identity 
in Equation (8) we obtain 
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detfi' 



dctO _ 

L ' 0<X N <---<X 1 <L-1 



det [/iAi-j+iW] det ^/t Ai+J v+i-i,j] 

■■■<Ai<L-l 
det [h x N + l-j—N — l+j+i 

{x}) det ( KAi+JV+ i-;,,) 

0<Ajv<---<Ai<L-1 

^ det (/ifcj-AT-i+ijx} j det ^Kfc.j) 



l<k!<---<k N <N+L-l 



det 



( iV+L-1 

(27) 



fc=i 

Hence det ft' is Casoratian in {x}. Next, we need to show that a Casoratian 
determinant is a solution of the bilinear identities of discrete KP, but this 
requires a number of definitions which we outline in the next section. 

4. Continuous KP, Miwa variables and discrete KP. 

A standard introduction to the continuous KP hierarchy is [21]. Miwa 
variables are discussed in detail in [22] where further references to their 
applications are provided. The discrete KP hierarchy was introduced in [9] , 
and further studied in [23] and [14]. In this note, we follow the treatment 
in [15]. 

Continuous KP is an infinite hierarchy of intcgrable partial differential 
equations generated in Hirota's bilinear form by expanding the bilinear 
identity 

I Ak_ e «t-t',fe) T ^_ e ( fc -i)j T {t + e (fc-i)} = (28) 
Jk=k oa 2vri 

where k G P 1 , the contour integral is around the point at infinity ki„fi n it y £ 

P\ {*} = {ti,t 2 ,t 3 ,---}, Z(t,k) = ESitifc*. e^- 1 ) = {l,^,^,-}, 

{t ± e(fc -1 )} = {ti ± \,t2 ± ^i*3 =t " ' ' }• The simplest KP equation in 
the hierarchy is 

[Df + :w\ - 4D1D3 j t • t = (29) 

where Di is the Hirota derivative with respect to tj. For the precise defini- 
tion of Di and that of the notation r • r, see [21]. 
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Continuous and discrete Miwa variables. In [23], Miwa introduced 
two infinite sets of variables, the continuous variables {x} = {xi,X2, ■ • ■}, 
and the discrete (and integer valued) variables {rn} = {mi,m2, • • • }. and 
showed that setting 



'i = X! mi T (30) 
»=i •* 

transforms T-functions of continuous KP to r-functions of a hierarchy of 
bilinear difference equations, namely discrete KP, studied in detail in [14]. 
These variables are now known as continuous and discrete Miwa variables, 
respectively, 

Multiplicities. From Equation (30), one can see that the discrete variables 
{to}, where 6 Z are multiplicities of the continuous variables {x}. In 
other words, m; > 1 is equivalent to saying that X{ occurs m, times in {x}, 
or that there are m, continuous variables that have the same value Xi . 

Discrete KP is an infinite hierarchy of integrable partial difference equa- 
tions in an infinite set of continuous Miwa variables {x}, where time evo- 
lution is obtained by changing the multiplicities {to} of these variables. In 
this note, we are interested in situation where the total number of contin- 
uous Miwa variables is finite, and the sum of all multiplicities is N. In 
this case, the discrete KP hierarchy can be written in bilinear form asnxn 
determinant equations 



' 1 x x 
1 x 2 

det 

k 1 x n 

where 3 < n < AT, and 



irVi{i}r.i{i} i 
x^ 2 t +2 {x}t^ 2 {x} 



2 T +n {x}T_ n {x} , 



= 



(31) 



'+* X^S — ' l J 'l i * * * t i i ' ' ' i TV J ' 

r_,{x} = r{x< mi+1) ,-- - ,x { r\--- ,4 m " +1) } (32) 

In other words, if t{x} has m; copies of the variable x^, then r+i{x} has nii+ 
1 copies of Xi and the multiplicities of all other variables remain the same, 
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while T~i{x} has one more copy of each variable except X{. Equivalcntly, 
one can use the simpler notation 



The simplest discrete KP bilinear difference equation, in the notation of 
Equation (33), is 

Xi(xj - x k )T{m,i + l,m j ,m k }T{m l ,m j + l,m k + 1} 
+ Xj(xk ~ x i )T{m l ,m j + l,mfc}r{m t + l,mj,mk + 1} 
+ x k (xi - Xj)T{mi,mj,m k + l}r{TOj + l,m,j + l,m k } = (34) 

where {xi,Xj,x k } G {x} and {rm,mj,m k } G {m} are any two (correspond- 
ing) triples in the sets of continuous and discrete (integral valued) Miwa 
variables. Equation (34) is the discrete analogue of Equation (29). 

Discrete time evolution in discrete KP. Each continuous Miwa vari- 
able Xi corresponds to a time variable in discrete KP. Time evolution in 
discrete KP, in direction Xi, is given by the discrete changes in the multi- 
plicities rrii of x t . Notice that as a multiplicity nii changes by ±1, the rank 
of the matrix Mj±i, where detMi±i = Tj±i remains the same as the rank 
of M, where det M = t. 

5. Bethe scalar products are discrete KP r-functions. 

In this section, we adapt the general treatment of [15] to the specific case of 
Slavnov's determinant expressions. We do this in detail to show explicitly 
that the general (and slightly abstract) identities and theorems in [15] apply 
to Slavnov's expressions. 

Re-arranging the elements of Slavnov's determinant. Given the 
NxN matrix f2' with elements 



T +i {x} = r{mi, • • • , (mi + 1), • • • , m N } } 
r_i{x} = r{(mi + !),-•• ,m;, • • • , (m N + 1)} 



(33) 



N+L-l 




(35) 



k=l 



let us consider the matrix 17" with elements 
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N+L-l 



^ Cik h k -j{x} 



(36) 



k=i 



which is obtained from Q' by reordering the rows of the latter from bottom 
to top, changing the rows and the columns and setting at — Hki- Notice 
that we use u) rather than lo" for the elements of f2" to simplify the notation. 
Since detfl' = (— 'j N ( N - 1 )/ 2 d e t fl 7 it is sufficient to show that Q satisfies 
the difference bilinear identities of discrete KP. 

Identities for the elements u>ij. It follows from Equations (4-5) that 
the elements w%j of Q" satisfy analogous identities 



u} ij {x 1 , 



.,X N } = Uij{xi, . . .,x N } 

+ XmUJij+^X!, . . . ,X$ , . . . ,X N } (37) 



(x r Xg) LU<ij {x\ , . • . , ^ , Xg \ . . . Xj\r} — 
,(2) 



x r Uij{xi, ...,x\.',.. .,x N } - x s Uij{xi, ...,x ( s ',.. .,x N } (38) 
From Equation (7), we see that 

A m u>ij{xi, . . -,x N } = Uij + i{xi, . . .,x N } (39) 

which is equivalent to the statement that det is Casoratian. 

Notation for column vectors with elements «y . We need the column 
vector 



u>i = 



f r ( m i) 
r ( TO i) 



(mi) 



' ) X N J 
■ ) X N 5 



..,X 



(m N ) 
N 



}) 



(40) 



and write 
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3 



( t (mi) ("»fc„+l) (tojv)-i \ 

,...,x ki ,...,x kn ,...,x N } 

,, r~("»i) J m ki+ 1 ) „(™*„+l) „(mjv)-i 

^2j\-ii 7---j x fci '•••' x fe„ f-^N I 

,,„.f>l) J" 1 *^ 1 ) „.("»*„ + 1 ) „,(™«h 

V u 'Wjl J 'l '•••i' , 'fc 1 '•••' X fc„ :---: x AT ]) 



(41) 



for the corresponding column vector where the multiplicities of the variables 
Xfcj , • • • , Xk n are increased by 1. 

Notation for determinants with elements U3ij. We also need the de- 
terminant 



r = det ^uj\ UJ2 • • • wjv j = | wi CJ2 • • • wtv | (42) 
and the notation 

r [fcl,...,fc»] = | ^.-.fcn] . . . ^.....fcn] | (43) 

for the determinant with shifted multiplicities. Next, and closely follow- 
ing [15], we derive two identities involving Casoratian determinants with 
elements Uij. 

Casoratian identity 1. The first identity that we need is 

X\- 2 T [1] = | Wi LJ 2 ■■■ LON-1 ^N-n+2 | ( 44 ) 

which is derived as follows. From Equation (43), we have 

Jl] I ,-j[l] -[l] ,-j[l] I ( A r\ 

In view of Equation (37), subtracting xi times column j + 1 from column 
j in this determinant for j = 1, 2, . . . , N — 1 allows us to write 

T [1] = | CJi UJ2 ■■■ W$ | (46) 

Multiplying column N by x\ and adding column N — 1 to the result, we 
obtain 



2:1 t [11 = | wi co 2 ■ ■■ wjv-i 4-i I (47) 
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Similarly, multiplying column N in Equation (47) by x\ and subtracting 
column N — 2 yields 

x\ t [1] = | coi w 2 ■■■ Wjv-1 I ( 48 ) 

Iterating this procedure by multiplying column N by xi and subtracting 
column N — j, we obtain Equation (44). 

Casoratian identity 2. The second identity that we need is 
II (x r - x s )r^-^ = 

l<r<s<n 

I - - ->["] ->["-!] -[1] I /^nl 

| Wl ... UJ N - n 0J l N L n+1 0J l N -n+l ■ ■ ■ ^N-n+1 \ I 49 ) 

which is derived as follows. From Equation (47), it follows that 

„ --[1,2] _ I rjP] ,-tP] ,-jI 2 ! ,-s[1.2] I 

X\ T — I w l w 2 w jV-l W N-1 | l tlU J 

which we can rewrite by subtracting x 2 times column j + 1 from column j 
for j = 1,2,..., N -2 as 

fl 21 I - - -> -[2] -.[1,2] I /,-, \ 

II T 1 ' 1 = | U)i W 2 ' ' ' Wat_ 2 W] V J _ 1 wV-i I ( 51 ) 

Multiplying column TV by (x\ — x 2 ) and applying Equation (38), we see 
that 

I \ fl 21 I - - -[2] -[1] I 

(Xi - X 2 )X 1 T l 1 = X 1 I Wi W 2 • • • CJAr_2 Wjv-1 W V-1 I 

- X 2 | Wl W 2 • • • WjV-2 ^N-l I ( 52 ) 

Since the last two columns of the latter determinant are identical, we obtain 

(.Tl - X 2 ) T [1 - 2] = j Wl ... WjV-2 U [ N-V <*>jV-l I ( 53 ) 

which establishes Equation (49) for n = 2. Now suppose inductively that 



II (*r - X.) T^-^ 



1 < r < s < n 



I -, -> . ; f ->[n — ll 1 .. - 1 , 

wi ... u N - n wy_ n+1 oj 1 n > ■ ■ ■ wy_ n+1 (54) 
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then analogously to Equation (47) , we have 

n 

1 1 < , n - x ») t[i "" m = 

n| -. -. _[n] -.[n-1] ^[1] 



2—1 l<r<s<n 



I _ _ _[n] ->[n-l] -.[1] I /rr\ 

wi ... u N - n wJ v J _„ w]y-„ ••• wV_ n (55) 



It follows that 



i—1 1 < r < s < ?i 

I ,-t[n+l] ,-!["+!] ,-![","+!] -.[n-l,n+l] ,-t[l,n+l] I . . 

I ->[n+l] -.[n,n+l] -»[n-l,n+l] -.[l,n+l] I / cc \ 

| wi ... cj jv _„_ 1 cj] v _„ j u l N _ n tv l N _ n ■■■ uj l N _ n | (56) 

Using Equation (38) repeatedly gives 



] [ {xi - x n+1 ) x 



Ki<r 



I _ _ ->[n+l] -.[n,n+l] -.[n-l,n+l] -.[l,n+l] I 

n 

nl -> - -»[n+l] -.[n] ->[n-l] -.[1] I /e^n 

i=l 

Combining this with Equation (56) shows that 



l<r<s<n+l 

|W1 ... WjV-n-1 W^!_„ w]v-„ W V-n| (58) 

thereby completing the proof of Equation (49). We are finally in a position 
to complete the proof that Slavnov's determinant expressions are discrete 
KP r-functions. 



Bilinear identities from Laplace expansions. Following [15], we con- 
sider the 2Nx2N determinant, which is identically zero, 
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(-W n n -[n] -[2] "\ 

UJi ■■■ WAr_i ^ N _ n+2 Ui • • • UjV-n+1 ^Ar_ Il+2 ' ' ' ^N-n+2 = q 
Oi • • • Ojv_i W™_ n+2 UJi ' • ■ Wat_„ +1 <3&L n+2 • • • J 

(59) 

where we have used subscripts to label the zero elements with the positions 
of the columns that they are in for notational clarity. Performing a Laplace 
expansion of the left hand side of Equation (59) in Nx N minors along the 
top NxN block, we obtain 

n 

i/=i 

| Wl • • • WjV-n+l^-n+2 • • • w) V -r l +2 a; W-r l +2 ' ' ' ^N-n+2 \ = ( 60 ) 

Using Equations (44-49), Equation (60) becomes 

n 

^(-)- 1 ^- 2 T M [] (^ - ^Jrl 1 --*--"] = (61) 

i/— 1 l<r<s<n 

which we recognise as the cofactor expansion of the determinant in Equation 
(31) using the last column. Hence we conclude that Slavnov's determinant 
expression for the XXZ Bethe scalar product is a r-function of discrete KP. 

6. Remarks. 

Shifted r-functions are not Bethe scalar products. A Bcthe scalar 
product that involves m, rapidities Xi, for i = 1,2, ••• ,i m ax, is a Ca- 
soratian determinant of a matrix of rank r = YL % i^T m «- Let us de- 
note the corresponding r-function by r — r{x^ mi \ • • • ,x\ mz \ ■ ■ ■ ,x^ N ^}. 
Now let us consider a time evolution of the latter, for example Tj+i = 
r{x^ ni \ ■ ■ ■ ,x| m,+1 \ • • • ,x^ N ^}. Time evolution has increased the mul- 
tiplicities by 1, but kept the rank of the corresponding Casoratian deter- 
minant the same, thus we cannot interpret rj+i as a Bethe scalar product 
and it remains unclear to us how to interpret the discrete time evolution of 
a Bcthe scalar product in the language of the XXZ spin chain. 
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Fermionization remains valid. Continuous KP r-functions can be writ- 
ten as expectation values of charged free fermion operators [21]. This re- 
mains the case for discrete KP r-functions and was the starting point of the 
results of [14,23]. In [13], the fermion expectation value version of Slavnov's 
determinant expression was obtained based on an earlier result [26]. It is 
straightforward to show that this result remains the same as the continuous 
KP time variables are restricted to be power sums of a finite and smaller 
number of continuous Miwa variables. 

Relation to the work of Krichever et al. As mentioned earlier, our 
result is close in spirit to that of Krichever et al. [7, 8] and works that 
followed including [24,25]. The starting point of [7] is that the Bethe eigen- 
values satisfy a bilinear identity that has the same structure as Hirota's 
bilinear difference equation and hence can be identified with r-functions of 
a discrete hierarchy. From this, a large number of interesting results follow, 
including an identification of the fusion rules of the transfer matrices of the 
quantum spin chain with Hirota's difference equations, that each step in 
the nested Bethe Ansatz approach to the spin chain is identified with a clas- 
sical Backlund transformation, and most interestingly that the eigenvalues 
of Baxter's Q operator are classical (suitably normalized) Baker- Akhiezer 
functions. On the other hand, our result is that it is the Bethe scalar prod- 
uct of a Bethe eigenstate rather than the corresponding Bethe eigenvalue 
that is identified with a discrete KP r-function, and we are far from ob- 
taining further results that are analogous to those of [7] . We hope that our 
identification is compatible with and complements that of [7]. 

Relation to the work of Sato and Sato. Equation (61) also follows 
from Theorem 3 of Sato and Sato [27]. We didn't know this when we 
obtained our proof, and the existence of more than one proof can only shed 
more light on the result obtained. 
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